Inert states of spin-S" systems 
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We present a simple but efficient geometrical method for determining the inert states of spin-S 1 
systems. It can be used if the system is described by a spin vector of a spin-S 1 particle and its energy 
is invariant in spin rotations and phase changes. Our method is applicable to an arbitrary S and it 
is based on the representation of a pure spin state of a spin-S particle in terms of 2S points on the 
surface of a sphere. We use this method to find candidates for some of the ground states of spinor 
Bose-Einstein condensates. 



States which are stationary points of the energy re- 
gardless of the exact form of the energy functional are 
called inert states A stationary point may change, 
say, from a local maximum to a global minimum as the 
parameters characterizing the energy vary, but for all pa- 
rameter values the inert state remains nevertheless a sta- 
tionary point of energy. The independence of the details 
of the energy functional is related to the fact that inert 
states are determined by the symmetry group of the en- 
ergy and the symmetry properties of the order-parameter 
alone. Inert states have been studied in the context of 
superfluid 3 He, where the analytical minimization of the 
energy functional is a very complicated task, but the in- 
ert states can be calculated straightforwardly [J, 0] , see 
also Q. It should be noted that sometimes the form of 
the ground state may depend on the parameters of the 
energy, in which case the ground state cannot be an in- 
ert state. It is therefore possible that not all stationary 
points of the energy are inert states, which is the case 
in 3 He [H, 0|. For spin-S systems in general, however, 
the calculation of inert states can also be tedious, as re- 
cently demonstrated for integer spin- values S = 1 — 4 [B[ , 
and the task seems to become increasingly complicated 
as the value of S increases. In this paper we present a 
simple geometrical description that allows one to deter- 
mine the inert states of spin-S systems. The method is 
applicable regardless of the value of spin and by using it 
we explicitly construct numerous inert states. Almost all 
inert states for S = 1,2,3 systems are ground states of 
the presently realized spinor Bose-Einstein condensates 

0, HSU la EES, which suggests that such states are 



also ground states of spinor condensates when S is larger 
than three. 

Our way of calculating inert states is based on a geo- 
metrical representation of spin states. It is well known 
that a pure spin state of a spin-I/2 particle can be char- 
acterized by a point on the surface of the Bloch sphere. 
Moreover, a pure spin state £ of a spin-S particle can 
be written in terms of 2S pure spin states of spin-1/2 
particles H, 13|. Combining these results, we see that 



£ can be expressed by 2S points on the surface of the 
Bloch sphere, which provides the geometrical description 
we are interested in. This decomposition was originally 
used to study the motion of a spin-S particle in a mag- 



netic field [I2j, |13j] . Recently it has been used to calculate 



the symmetries of the ground states of spinor condensates 

0,13, 



16, V, 



with S = 2 and S = 3 [14|, UM, see also Refs - 
This decomposition is well suited for such purpose, since 
it shows explicitly how spin vectors change in rotations. 
Rotating £ corresponds to rotating the point configura- 
tion describing this vector with respect to the same axis 
and through the same angle. Let £ = ^] M= _ g £m|S, M) 
be a spin vector of a spin-S particle. The points on the 
surface of the Bloch sphere characterizing £ can be ob- 

HQ 



tained from the equation 
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In the second expression we have chosen A € C in such 
a way that we can define a k = cos(#fc/2)e~ 4Vfc / 2 and 
f3 k — s\\Y{9k/2)e lipk / 2 for all k. Every pair (ctk,Pk) deter- 
mines a point on the Bloch sphere located at (9k, <fik)- In 
the rest of the paper we call the polynomial of Eq. JT]) the 
characteristic polynomial of £. If we are given two sets of 
points, then the characteristic polynomial of the configu- 
ration consisting of both sets is obtained by multiplying 
the characteristic polynomials of the two sets. Equa- 
tion ([T]) shows that the point distribution is the same 
for £ and ht;, where h is an arbitrary non-zero complex 
number. 

We now move on to calculating the inert states. We 
assume that the system is described by a spin vector 
of a spin-S particle, denoted by £, and that the energy 
functional is invariant in rotations and phase transfor- 
mations of £. Thus the symmetry group of the energy 
is G = U(l) x SO(3), where SO(3) acts via its irre- 
ducible 2S + 1 dimensional representation and the ac- 
tion of £7(1) corresponds to multiplication by a com- 
plex number of modulus one. We define the isotropy 
group of £ by = {g e G \ g ■ £ = £}, i.e., it consists 
of those elements of U(l) x SO(3) which leave £ fixed. 
Two isotropy groups, say and H^, are conjugate if 
= gHsjg^ 1 = {ghg^ 1 \ h S H^} for some g e G. 
In Ref. [HI it has been shown that £ is an inert state 
if, for every £' infinitesimally close to £, and H^i are 
not conjugate. Here it is assumed that £' cannot be ob- 
tained from £ by a rotation and a phase change and that 
1 1 £11 = IICII- In the present case each £ determines some 
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point configuration on the surface of a sphere. A small 
variation in the positions of the points leads to a small 
change in £. Therefore, if all small changes necessarily 
lead to a new isotropy group which is not conjugate to 
i?f , then £ is an inert state. The elements of the isotropy 
group consists of pairs (u,r), where u 6 U(l) and 
r € 50(3). The latter are the elements of some sub- 
group K of 50(3). K can be determined by finding out 
all rotations which leave the point configuration invari- 
ant. We call K the symmetry group or symmetry of the 
point configuration. It should not be mistaken for the 
isotropy group, see Ref. [l7j |. 

We find the inert states by determining point configu- 
rations which cannot be modified without changing the 
symmetry group K. However, even if K is invariant 
as the positions of the points are varied, it is possible 
that the elements u £ U(l), and therefore the isotropy 
group, change. Let s be some symmetry axis of the point 
configuration. Because of the rotational degree of free- 



dom, we can choose s = z 
for some a rot - Writing £ = ^ 
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£m\S,M) we get 
' £m for all M. Assume that p n (p s ) of 
the points characterizing £ are at the north (south) pole. 
Then Eq. fl} shows that the highest order term in the 
characteristic polynomial is x 2S ~ p " and the lowest order 
term is x Pn . Therefore £s-p s , £-S+p„ an d we get 



Pn)otrot 

- S)a rot 



(mod 2ir) 
(mod 2-71"). 
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A similar result has been presented in Refs 
(u, r) S -ffj and the rotation angle of r is a ro t 



(2) 

3- If 

then 



u = e . We see that if the elements of the 50(3) part 
of remain unchanged as the points move, the U(l) 
parts can change only if the number of points on a sym- 
metry axis changes. If the smallest non-zero rotation 
angle leaving the points fixed is a ro t — lixjn, then the 
only way to move points from the symmetry axis without 
affecting the symmetry group is to remove symmetrically 
kn points from the axis. Here k = 1, 2, . . .. Equation @ 
shows that in this case the U(l) parts of remain un- 
changed, so it is possible to change the phase terms only 
if also the symmetry group changes. When calculating 
inert states we can thus concentrate only on the changes 
in the symmetry group. 

The possible symmetry groups are the continuous 
groups 0(2), 50(2), and the finite groups /, C n ,D n , T, O 
and Y. I is the trivial group consisting of the identity 
element alone, C n (D n ) is the cyclic (dihedral) group of 
order n, T is the symmetry group of the tetrahedron, O 
that of the octahedron or the cube, and Y that of the 
icosahedron and the dodecahedron. We now calculate 
which of these subgroups determine inert states. We first 
find point configurations that have some given group as 
their symmetry group, and then see which of these can- 
not be modified without changing the symmetry group. 



In this context Ref. [19| is very useful. In the following 
we choose the z-axis to be one of the symmetry axes. 

Orthogonal group 0(2). The point configuration has 
to be invariant in arbitrary rotations about the z-axis 
and in a rotation through ir about an axis in the xy- 
plane. Therefore there has to be an equal number of 
points at the north and south poles and no points else- 
where. Clearly any change in the relative positions of the 
points will change the symmetry group, so this state is 
an inert state. Now (8k, fk) = (0,0) for k = 1,...,5, 
(8k, fk) = (tt, 0) for k = S + 1, . . . , 25, and the spin vec- 
tor obtained from Eq. (p]) is £0(2) = 1^0). This vector is 
clearly possible only if 5 is an integer. Here and in the 
rest of the paper the subscript of £ gives the symmetry 
group of the point distribution and the superscript gives 
the value of the spin. 

Special orthogonal group 50(2). This symmetry is 
present if there is only one symmetry axis and the point 
configuration is invariant in arbitrary rotations about 
this axis. Thus all the points have to be at the poles. 
The number of points at different poles cannot be equal, 
because then the configuration would be invariant un- 
der 0(2). Also now the symmetry group changes if the 
relative positions of the points are changed, so this is 
an inert state. If there are 5 + M points at the north 
pole, then (8 k ,ip k ) = (0,0) for k = 1,...,S + M and 
(8 k , <p k ) = (tt, 0) for k = S + M + 1, . . . , 25. These give 
£fo(2) = I^-^Oj where M > 0. It is sufficient to con- 
sider positive M only, since |5, M) and |5, — M) can be 
obtained from each other by a rotation. Now inert states 
can exist also when 5 is a half integer. 

Cyclic group C n . The only symmetry axis is the z- 
axis, and the system is invariant in rotations through 
27rfc/n about this axis. There has to be a point at 
(8, ip) = (8c, <fc + 2nk/n) for some fixed 8c S (0, ir), 
ipc £ [0, 2ir/n), and for each k = 1, ...,n. There may 
be points also elsewhere, as long as the configuration 
is invariant under C n . The value of 8c can always be 
changed by some 88c without altering the symmetry 
group. This change cannot be achieved by a rotation. 
Therefore the cyclic groups do not determine any inert 
states. In Fig. [lja) we show an example of the O5 sym- 
metry in an 5 = 5/2 system. 

Dihedral group D n . Because C n is a subgroup of D n , 
there has to be a point at (8, ip) = (8b, Pd + 2nk/n). 
Now there are n symmetry axes in the xy-plane, so there 
is an equal number of points at the north and the south 
pole. If there are points outside the xy-plane and the 
poles, it is possible to change the ^-coordinate of these 
points without changing the symmetry group, so these 
configurations do not give inert states. In Fig. |Tf b) 
this is illustrated for a configuration with D4 symme- 
try. The case where there are points only in the xy- 
plane and at the poles can be divided into three cases. 
Below k = 1,2, ...,n, / = 1,2,..., and ipi S [0, 2ir/n) 
(I.) n or more points at each pole and In points in the 
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xy-plane. Now n points can be moved from each pole 
towards the equatorial plane without changing the sym- 
metry, see Fig. QJa). (II.) Less than n points at each 
pole, two or more points at (9, ip) = (ir/2, ipi + 2itk/n), 
and In points elsewhere in the xy-plane. It is possible to 
move one point up and one down by the same amount 
at each (tt/2, ip\ + 2 / nk/n) without changing the symme- 
try group. (HI.) Less than n points at each pole and 
one point at each {9, ip) = (ir/2,ipi + 2itk/n), ipi ^ ipy 
if I ^ When the symmetry axes in the xy-plane are 
taken into attention, it turns out that / = 1 is the only 
possibility (20j . 

We therefore conclude that an inert state with D n sym- 
metry is possible only if there is one point at (9, ip) — 
(7r/2, 2nk/n) for k — 1,2, ...,n and p < n points at 
each pole (we have chosen ipi =0). It follows that in 
a spin-5 system there is an inert state corresponding to 
D n only if (25 + 2)/3 < n < 2S. Furthermore, the con- 
dition n = 2S — 21, where / is a non-negative integer, 
has to hold. This is because, for fixed 5, n can be de- 
creased only by putting an equal number of points at each 
pole. Therefore for integer (half integer) spin only even 
(odd) values of n are possible, the smallest value of n 
being given by the smallest even (odd) integer bounding 
(25 + 2)/3 above. It is to be noted that if 5 = 3, n = 4, 
the symmetry of the point configuration is not -D4, but 
that of an octahedron. If 5 = 2, n — 2 or n = 4, the 
symmetry is -D4. The n = 4 case corresponds to a square 
in the horizontal plane, while n = 2 gives a square in a 
vertical plane. We define ?£ 2S _ 21 = \S,S-l) + \S,-S+l). 
Using Eq. JT]) it is easy to verify that this is a vector with 
D2S-21 symmetry. Dihedral inert states exist if 5 > 1 for 
both integer and half integer 5. 




FIG. 1: Some non-inert states. In (a) 8 and in (b) the height 
of the box can be changed without affecting the symmetry. 
In (c) we show the truncation of an octahedron. 




FIG. 2: Examples of inert states. In (a) cube and octahedron 
and in (b) cube and tetrahedron inscribed in a sphere simul- 
taneously. In (a) the symmetry is that of octahedron/cube 
and in (b) that of tetrahedron. In (c) icosahedron and in (d) 
dodecahedron. 

The remaining possible symmetry groups are tetrahe- 



dral, octahedral and icosahedral groups. We calculate 
only those inert states which are obtained by placing 
points at the vertices of the convex regular polyhcdra. 
The latter are often called the Platonic solids, and they 
are the tetrahedron, octahedron, cube, icosahedron and 
dodecahedron. Any change in the relative positions of 
the vertices changes the symmetry group (assuming that 
rotations are excluded). It is possible that there are also 
other ways to distribute points on the surface of a sphere 
with, say, octahedral symmetry than to place them at the 
vertices of the octahedron or cube. However, these distri- 
butions correspond to polyhedra which are not regular, 
and therefore the points can probably be moved without 
changing the symmetry group. Platonic solids have an 
even number of vertices, which means that inert states 
with half integer 5 are not possible. 

Octahedral group O. Octahedron and cube can 
be inscribed in a sphere simultaneously so that the 
combined point configuration has octahedral symme- 
try, see Fig. \T[ a). This is obtained by putting the 
vertices of the cube at (9,ip) = (9c u be,^k/2), (n — 
9c'ube,'xk/2), and those of the octahedron at (9,p) — 
(0,0),(7r,0),(7r/2,7r(l + 2fc)/4). Here k = 1,2,3,4 and 
9cube = tan _1 (V2)- The corresponding characteristic 
polynomials are x 8 — 14a; 4 + 1 and x 5 + x for the cube 
and octahedron, respectively. The vectors obtained from 
these polynomials are & ube = V5|4,4) - \/l4|4,0) + 
\/5|4, -4) and & cta = |3, 2) + |3, -2). If there are m (n) 
points at each vertex of a cube (octahedron) , the charac- 
teristic polynomial becomes (a; 8 — 14x 4 + l) m (x 5 + x) n . 
Corresponding spin vector can be calculated by compar- 
ing this polynomial with the first expression in Eq. fT]). 
For example, if m = n = 1, the characteristic polynomial 
becomes x 13 — 13x 9 — 13x 5 +x, which gives ^cube+Octa = 
\/lT]Y, 6) - VT3|7, 2) - Vl3|7, -2) + \ZTT|7, -6). It is im- 
portant to notice that if there are four or more points at 
each vertex of the octahedron, it is possible to truncate it 
continuously, see Fig. QJc). Truncation does not change 
the symmetry or isotropy group but it changes the state, 
so the octahedron determines an inert state only if there 
are less than four points at each vertex. Similarly the 
cube can be truncated if there are three or more points 
at each vertex. Combining these results, we see that in- 
ert states determined by the cube and octahedron are 
possible if 5 = 3m + An, where m + n > 1, m < 3, and 
n < 2. 

Tetrahedral group T. Tetrahedron has four 
vertices, which we choose to be at (9, cp) = 

{9Tetra , 7r/4), (9 T etra , 5"7r/4) , (n - 9 T etra , 37r/ 4), (% - 

&Tetra,7Tr/4:), where 9xetra = tan" - l (V2). The char- 
acteristic polynomial is x 4 + 2i^/3x 2 + 1, which gives 
= 1 2, 2) + iV2\2,Q) + |2,-2). The tetrahedron, 
octahedron, and cube can be inscribed in a sphere 
simultaneously, see Figs.EJa) and[2]Jb). The vertices can 
be chosen to be at the points given above. The resulting 
configuration has tetrahedral symmetry, which therefore 
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TABLE I: The inert states for 5" = 1 - 4. The only symmetry 
group that is not manifested here is the icosahedral group. 



Spin 


Inert states 


S = 1 


£ SO (2) = |M>: e O (2) = |i,0) 


5 = 2 


£ so(2) = |2, 2} and j2, 1), £d 4 = \2, 2} + |2, -2), 
?0(2) = |2, 0), and fret™ = |2, 2} + i-y/2|2, 0} + |2, -2} 


S = 3 


Cso( 2) = |3, 3), |3, 2), and |3, 1), f 0(2) = |3, 0), 
£d 6 = |3, 3) + |3, -3), and £ 0c ta = |3, 2) + |3, -2} 


S = 4 


£so(2) = |4, 4), |4, 3), |4, 2), and |4, 1), £ 0(2) = |4, 0), 

fr 8 = |4, 4) + |4, -4), £ D „ = |4, 3) + |4, -3) 
£d 4 = |4, 2) + |4, -2), fr ctra = V7J4, 4) + 2iV3|4, 2} 
— v/I5|4, 0) + 2zV3|4, -2) + V7|4, -4). 
£c«i,e = \/5|4, 4} - x/T4|4, 0) + v/5|4, -4) 



can be present at least if 5 = 21 + 3m + 4n. Here I > 1 
and m,n > 0. Because of the restrictions imposed by 
truncation, the conditions I + n < 2 and m < 3 have to 
hold. 

Icosahedral group Y. Icosahedron can be truncated 
if there are five or more points at each vertex, while 
for dodecahedron this is possible if there are more than 
two points at each vertex, see Figs. [2jc) and [2fd). 
Thus these symmetries are possible if 5 = 6m + lOn, 
where m + n > 1, m < 4, and n < 2. We obtain 
e ico = %/7|6,5) - VTT|6,0) - V7|6,-5) and ^° ode = 

17|10, 10) + v/57|10, 5) + y^|10, 0) - %/57|10, -5) + 

17| 10, —10), details of the calculation will be presented 
in Ref. 

In Table Q] we list the inert states for integer-valued 
spin 5=1-4. All except fo 4 = |4, 2) + |4, -2) have also 
been obtained using a different method [5]. This state 
was omitted in Ref. [f| because it was assumed that the 
isotrop y gr oup of £d 4 is a subgroup of the isotropy group 
°f £d 8 [21(. But we find that this is not the case, as can 
be seen by considering a rotation through it/2 about the 
z-axis. In this case Eq. ^ shows that for D4 symmetry 
6 = 7T, while for D% symmetry 6 = 0. 

It is interesting to compare the inert states of Table Q] 
with the ground states of spinor Bose-Einstein conden- 
sates with 5 = 1,2,3. For 5 = 1 there are no other 
ground states than the inert states [1, 0|- If 5 = 2 
it seemed initially that also non-inert ground states are 
possible due to an accidental degeneracy [H, [2§| , but re- 
cently it has been shown that due to quantum and ther- 
mal fluctuations this degeneracy is lifted, after which 
only inert ground states can occur [2^, [2J]. However, 
there is one inert state, namely |2,1), which is never 
a ground state. If 5 = 3, the inert states are ground 
states but also numerous other ground states are possi- 
ble @,[I3,Ii3- These results suggest that the inert states 
for 5 > 3 are ground states of spinor condensates at 
least for some values of the scattering lengths. This has 



interesting consequences considering topological defects, 
as it shows that monopoles made possible by the |5, 0) 
ground state and non-Abelian vortices arising from the 
-D„-symmetric states are possible in spinor condensates 
regardless of the value of spin. In this paper we have not 
considered the time-reversal symmetry. Its effects will be 
discussed elsewhere [20| . To conclude, we have presented 
a simple geometrical method to calculate the inert states 
of systems which are described by a single particle wave 
function of a spin-5 particle and which have U(l) x 50(3) 
invariant energy functional. The possible inert states can 
be identified and classified by their symmetry groups. 
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